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ADSORPTION,  DIFFUSION,  AND  EPITAXY  OF  SI 
ON  SI(IOO)* 


ZHENYU  ZHANG  and  HORIA  METIU 


Department  of  Chemistry  and  Physics 

University  of  California 

Santa  Barbara,  CA  93106,  U.S.A- 

ABSTRACT;  We  review  oxir  recent  work  on  the  migration  and  the 
self-organization  of  the  Si  atoms  deposited  on  a  Si(100}-(2X1)  surface. 
The  silicon  silicon  interactions  are  modelled  by  the  Stillinger-Weber 
potential,  and  we  use  molecular  dynamics  and  Monte  Carlo 
simulations,  total  energy  calculations,  as  well  as  phenomenological 
approaches  to  study  the  following  topics  of  current  interest:  the 
effects  of  surface  reconstruction-unreconstruction  on  the  adsorption 
dynamics;  the  spatial  distribution  of  the  adsorbate  in  the  early  stages 
of  deposition;  the  migration  of  an  adatom  on  a  terrace  or  across  a 
step;  the  role  of  adatom-adatom  interactions  in  the  formation  of 
anisotropic  islands;  the  transformation  of  a  surface  with  single  layer 
steps  to  one  with  double  layer  steps;  the  effects  of  kink-kink 
interaction  in  step  roughening,  etc.  Wherever  possible,  the  results 
obtained  are  compared  with  the  findings  of  recent  experiments;  we 
find  that  our  qualitative  conclusions  are  in  agreements  with  the  STM 
observations. 
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L  INTRODUCTION 


For  technological  and  scientific  reasons  one  would  like  to  have  a 
better  understanding  of  the  kinetic  processes  taking  place  during 
epitaxial  deposition  on  semiconductor  surfaces.  The  bulk  of  our 
information  comes  from  electron  scattering  measurements^'^ 
performed  in  situ  during  growth.  If  the  deposition  is  made  on  a  fiat 
surface  the  main  issue  is  whether  the  growth  is  two  or  three 
dimensional.  Two  dimensional  growth  means  that  the  crystal  is 
grown  layer  by  layer.  During  the  deposition  one'atom-high  islands 
grow  and  join  each  other  to  form  a  surface  layer,  before  a  significant 
number  of  atoms  will  be  trapped  on  top  of  them.  Three  dimensional 
growth  is  opposite.  New  one-atom-high  islands  grow  on  older  one- 
atom-high  islands,  which  have  grown  on  older  islands,  etc.  Four  to 
five  layers  are  simultaneously  being  completed.  Electron  scattering 
can  distinguish  between  these  two  growth  regimes  because  they  lead 
to  a  different  evolution  of  surface  roughness.  During  three 
dimensional  growth  the  surface  is  always  rough;  during  two- 
dimensional  growth  the  roughness  of  the  surface  oscillates:  it  is 
highest  when  a  half  layer  is  deposited  and  lowest  whenever  a  new 
layer  is  completed.  Since  rough  surfaces  diminish  specular  (or 
diffractive)  scattering  (electron^’^  or  helium^),  the  specular  (or  any 
diffracted)  beam  intensity  oscillates  if  the  growth  is  two  dimensional. 

Ihe  substrate  surface  is  not  always  fiat.  Often  the  presence  of 
steps  is  unavoidable,  and  recently  vicinally  cut  substrates  have  been 
deliberately  used  for  the  growth  of  new  structures.  The  growth  on 
such  stepped  sxirfaces  has  new  features.  If  the  atoms  falling  on  the 
terrace  have  high  mobility  they  reach  a  nearby  step  before  they  have  a 
chance  to  meet  each  other  on  the  terrace.  If  they  stick  to  the  step,  the 
deposition  causes  step  propagation.^''^  In  the  opposite  limit  the  atoms 
form  islands  on  the  terraces.^'®  Step  propagation  causes  no  change 
in  surface  smoothness  or  in  the  periodicity  of  the  steps.  Therefore, 
the  diffraction  spots  due  to  terrace  periodicity  persist  throughout  the 
deposition  and  the  intensity  of  the  specularly  scattered  beam  is 
unaffected.  Two  dimensional  island  growth  on  terraces  causes 
atomic  level  roughness  on  the  terrace  until  the  deposition  of  a 


monolayer  is  completed.  The  diffraction  spots  characteristic  of  step 
periodicity,  and  the  high  intensity  in  the  scattered  beam 
(characteristic  of  surface  flatness),  deteriorate  as  the  coverage  is 

fi  n 

increased  and  recovers  when  one  layer  is  completed.  ’ 

In  the  case  of  multicomponent  systems  (e.g.  Ga  and  A1 
deposition  on  stepped  GaAs  surfaces)  the  step  can  induce  chemical 
segregation  (e.g.  A1  segregates  near  the  step  and  Ga  on  the 
remainder  of  the  terrace).®  One  would  like  to  know  when  and  how 
such  segregation  takes  place,  and  how  smooth  is  the  border  between 
the  phases.® 

From  the  point  of  view  of  the  atomic  theory  the  information 
provided  by  these  experiments  is  rather  crude.  Therefore,  it  is 
perhaps  most  reasonable  to  try  to  understand  them  by  using  kinetic 
simulations  in  which  the  adsorbed  atoms  tindergo  jumps  on  a  square 
lattice  and  interact  with  each  other  by  simple  rules.®*  The 
influence  of  the  finer  details,  such  as  the  reconstruction  of  the 
surface,  the  presence  of  multiple  adsorption  sites,  and  the 
complicated  nature  of  the  interatomic  potentials,  are  left  out  of  the 
investigation.  While  important,  such  simulations  are  not  reviewed 

here.^® 

The  development  of  the  scanning  tunneling  microscope  (STM) 
made  it  possible  to  investigate  surface  growth  at  an  atomic  level. 
STM  is  not  yet  an  universal  tool:  no  measmements  have  been  as  yet 
performed  under  the  high  temperature  conditions  used  in  epitaxy 
and  the  pictures  taken  on  some  of  the  practically  interesting 
materials  do  not  show  much  detail.  There  is  reason  to  believe  that 
these  shortcomings  are  temporary.  On  the  other  hand,  for  the 
surfaces  and  the  conditions  for  which  STM  works  well  the  quality  of 
the  information  is  astonishing. 

An  example  of  what  is  possible  is  provided  by  the  recent  STM 
studies®*^. 13, 16-22  of  the  Si(lOO)  surface  covered  with  adsorbed  Si 
atoms.  These  measurements  have  determined  the  adsorption  sites 
on  the  surface;  studied  the  shape  of  the  islands  formed  by  a  small 
number  of  atoms;  investigated  how  the  shape  of  the  steps  changes 
under  the  influence  of  heating  or  atom  deposition;  determined  the 
configuration  of  the  three  dimensional  islands;  and  provided 
complete  statistical  information  on  step  roughness. 


r 


Theory  can  explain  these  observations  only  by  using  models 
which  incorporate  the  atomic  details:  surface  reconstruction,  the 
location  of  the  adsorption  sites,  the  barriers  between  them,  realistic 
adsorbate  -  adsorbate  interactions,  etc. 

In  this  article  we  review  our  group's  efiForts23*28  to  provide  a 
qualitative  (or  semi-quantitative)  understanding  of  the  results 
obtained  by  STM  measurements  on  the  Si(lOO)  surface.  For  the  sake 
of  simplicity  we  deal  with  processes  that  either  involve  isolated 
adatoms,  or  low  coverage.  This  is  a  first  step  towards  a 
comprehensive  study  of  growth  kinetics. 

Because  of  the  length  of  the  article  it  is  usehd  to  siimmarize  our 
main  findings  and  explain  how  the  material  is  organized.  All  the 
results  presented  here  are  obtained  by  using  molecular  dynamics, 
Monte  Carlo  simulations  or  total  energy  calculations  with  the  silicon- 
silicon  interaction  proposed  by  Stillinger  and  Weber  (SW).29  As  will 
be  discussed  in  more  details,  there  are  many  other  empirical 
potentials  available  for  Si,  and  density  functional  and  tight  binding 
calculations  are  also  becoming  practical.  Unfortunately,  there  is  no 
proof  that  any  of  them  describes  all  silicon  properties  adequately.  A 
specific  version  obtained  by  fitting  a  specific  set  of  data  turns  out  to  be 
particularly  good  at  describing  several  properties  of  the  system,  but 
fail  with  others.  We  use  the  SW  potential  because  it  is  qualitatively 
reasonable:  it  reproduces  well  most  of  the  bulk  properties  of  solid  and 
liquid  Si,  and  many  of  the  properties  of  the  Si(lOO)  surface.  Our 
qualitative  conclusions  that  can  be  tested  experimentally  are  in 
agreement  with  the  observations.  Additional  information  on  the 
potentials  is  given  in  Section  2. 

In  Section  3  we  study  the  dynamics  of  single  atom  adsorption 
and  determine  the  absorption  sites.23,24  were  surprised  to  find 
that  there  are  all  together  thirteen  possible  absorption  sites,  within 
the  (1X2)  unit  cell  of  the  reconstructed  Si(lOO)  surface  shown  in  Fig.  1. 
Similar  calculations  based  on  the  density  fimctional  method30,31  or 
the  Tersoff  potential32,33  also  suggest  a  large  number  of  adsorption 
sites,  but  they  differ  in  the  precise  location  and  binding  energy.  At 
the  present  time  it  is  not  possible  to  conclusively  decide  which  gives  a 
better  description  of  the  real  system. 

One  might  naively  think  that  the  incoming  atoms  prefer  to  bind 


to  the  site  with  the  highest  binding  energy,  which  is  favored  by 
thermodynamics.  This  is  not  the  case.  In  our  simulations  we 
observe  a  highly  inverted  population:  the  site  with  the  highest 
binding  energy  is  not  reached  by  direct  absorption  and  one  of  the 
thermally  metastable  sites  (with  low  binding  energy)  has  the  highest 
occupation.  While  interesting,  this  phenomenon  does  not  affect 
growth  kinetics  since  site  to  site  hopping  restores  the  thermal 
population  within  microseconds. 

In  Section  4  we  estimate  the  site  to  site  hopping  rates  and 
discuss  the  nature  of  single  atom  migration  on  a  lattice  with  multiple 
sites  in  the  unit  cell.24,25  We  find  that  the  migration  of  an  atom  on 
the  Si(100)-2X1  surface  is  highly  anisotropic.  Qualitatively  it  can  be 
described  as  an  one-dimensional*  random  walk  on  top  of  a  substrate's 
dimer  row  (See.  Fig.  1),  with  brief  side  trips  to  a  site  located  at  the 
fiank  of  the  row. 

After  himdreds  of  thousands  of  such  random  jumps  the  atom 
can  get  trapped  at  a  site  with  high  binding  energy,  located  between 
the  rows.  This  migration  can  be  described  by  a  diffusion  equation 
only  on  a  time  scale  shorter  than  the  time  in  which  the  atom  is 
trapped.  The  diffusion  coefficient  is  a  complicated  function  of  the 
jumping  rates  between  the  many  sites  involved.  Since  each  rate 
constant  satisfies  an  Arrhenius  formula  with  a  different  barrier  and 
pre-exponential,  the  temperature  dependence  of  the  diffusion 
coefficient  does  not  have  an  Arrhenius  form.  However,  it  turns  out 
that  within  the  narrow  temperature  range  employed  in  crystal 
growth,  the  temperature  dependence  of  the  diffusion  coefficient  is 
well  fitted  by  an  Arrhenius  formula.  But,  except  for  special  cases, 
the  activation  energy  and  the  pre-exponential  differ  from  those  of  the 
site  to  site  jumping  rates. 

The  anisotropy  of  migration  predicted  by  our  calculations  is 
consistent  with  Lagally's  data.5.i3  and  Metiu^^fb)  presented 
simple  chemical  arguments  that  suggested  that  such  anisotropy 
exists  on  the  (100)  faces  of  all  materials  with  diamond  or  zincbiende 
structure. 

In  most  surface  science  systems  the  adsorbed  atoms  interact 
more  strongly  with  the  substrate  than  with  each  other.  Because  of 
this  the  hopping  rates  are  not  radically  altered  by  coverage.  This  is 


not  the  case  for  silicon.  The  studies  presented  in  Section  5  show  that 
the  rate  of  jumping  from  one  site  to  another  is  very  strongly  affected 
by  the  occupation  of  the  neighboring  sites.^®  This  dependence  is  not 
simple.  For  example,  it  is  not  a  fimction  only  of  the  number  of  atoms 
surroimding  a  given  site  but  it  strongly  depends  on  the  precise 
position  of  the  occupants. 


Fig.  1. 


The  STM  pictures  show  that  a  small  number  of  atoms  deposited 
on  the  surface  aggregate  into  thin  long  islands  perpendictilar  to  the 
substrate's  dimer  rows.l®'^'^  The  borders  of  these  islands  are 
unusually  smooth.  Moreover,  it  is  extremely  puzzling  to  see  that  the 


long  coast  of  the  island  is  perpendicular  to  the  direction  of  rapid 
migration.  This  means  that  even  though  the  single  atoms  on  the 
surface  are  most  likely  to  make  contact  with  the  long  coast  of  an 
island,  there  is  a  mechanism  that  transports  them  on  the  short  shore 
where  they  stick  and  increase  the  island's  length.  The  calculations 
reported  in  Section  5  find  that  an  atom  arriving  at  the  long  coast  is 
likely  to  dimb  up  on  the  island  and  travel  along  the  shore  towards  the 
island's  end.  In  this  trip  the  atom  can  either  fill  the  vacancies  along 
the  border,  making  it  smooth,  or  reach  the  end,  descend  onto  the 
substrate  and  stick  to  the  end  of  the  island  making  it  longer. 

The  motion  of  adsorbed  atoms  on  stepped  surfaces  are  of  great 
interest  because  they  are  often  used  as  a  substrate  in  epitaxy  (e.g.  to 

grow  tilted  superlattices®).  Often  steps  are  present  on  a  surface  even 
when  they  are  unwanted;  it  seems  that  it  may  be  impossible  to  cut  a 
Si(lOO)  surface  without  steps.  Finally,  even  if  the  substrate  stirface  is 
step  free,  subsequent  deposition  creates  islands  and  the  borders  of 
these  islands  are  steps.  For  these  reasons  we  have  studied  in  some 

detail  in  Section  6  the  motion  of  the  Si  atoms  along  and  across  steps.^^ 

Our  main  purpose  is  to  understand  a  striking  experiment 
performed  by  Hoeven  et  al^  They  deposit  half  a  monolayer  of  silicon 
on  a  Si(lOO)  surface  which  has  single-atom-high  steps.  The  absorbed 
atoms  migrate  to  form  a  surface  having  exclusively  double-atom-high 
steps.  To  do  this,  the  atoms  must  leave  every  other  terrace  and 
congregate  on  the  terraces  in  between.  Our  calculations  suggest  a 
kinetic  mechanism  for  this  process.  On  a  surface  with  steps  having 
a  single  atom  height  the  terraces  are  not  equivalent:®^  the  dimer 
rows  on  one  type  of  terrace  are  perpendicular  to  the  step;  on  the  other 
they  are  parallel  to  it.  These  terraces  alternate:  a  terrace  of  a  given 
type  is  always  sandwiched  between  two  terraces  of  the  other  type. 
The  atoms  landing  on  a  terrace  can  rapidly  reach  the  border  if  the 
dimer  rows  on  the  terrace  are  perpendicxilar  to  the  steps.  This 
happens  because  the  mobility  of  the  adsorbed  atoms  is  highly 
anisotropic  and  takes  place  along  the  dimer  rows.®'^3, 24*26, 30-32  Qur 

calcxilations®^  show  that  when  such  atoms  reach  the  border  there  is  a 
significant  probability  that  they  will  cross  it  and  move  on  the 
neighboring  terrace.  The  atoms  on  terraces  having  dimer  rows 
parallel  to  the  steps  bordering  them  tend  to  move  parallel  to  the  steps 


and  have  a  low  probability  of  reaching  the  border.  Therefore  they 
tend  to  remain  on  the  terrace  on  which  they  landed.  As  we  show  in 
Section  6  these  two  observations  explain  why  the  deposition  of  half  a 
monolayer  of  silicon  on  the  stepped  surface  turns  all  the  single  atom 
high  steps  into  double  atom  high  ones. 

The  last  Section  is  concerned  with  the  observation  that  the 
single-atom-high  steps  that  are  perpendicular  to  the  dimer  rows  on 
the  terrace  above  them  are  very  rough  and  have  a  large  number  of 
irregular  kinks.^®  Lagally’s  group^^  has  used  the  STM  pictures  of 
these  steps  to  determine  the  probability  of  finding  kinks  of  a  given 
height  and  a  given  distance  between  the  neighboring  kinks.  They 
have  also  shown  that  the  observed  distribution  can  be  fitted 
reasonably  well  by  a  model  which  assumes  that  the  energy  to  create  a 
kink  is  proportional  to  the  kink  length  and  is  independent  of  the 
nature  of  the  neighboring  kinks.  Here  we  show  that  the  agreement 
with  the  data  is  improved  when  we  take  into  account 

phenomenologically  the  kink  -  kink  interactions.^®  Moreover,  at  a 
more  microscopic  level  we  find  that  when  the  dimer  row  on  a  terrace 
are  perpendicular  to  the  step  below,  the  step  growth  is  rather 
peculiar.  To  a  first  approximation  each  dimer  row  grows 
independently.  This  lack  of  correlation  between  the  row  growth  leads 

to  step  roughness.^^ 

2.  MODEL  POTENTIAL  AND  METHODOLOGY 

Throughout  this  paper  we  use  the  potential  proposed  by 
Stillinger  and  Weber  (SW).^^  This  has  a  pairwise  part  whose  main 
function  is  to  describe  the  energy  needed  to  stretch  a  bond,  and  a 
three  body  interaction  whose  role  is  to  keep  the  bonds  at  the  correct 
angle.  Furthermore,  each  Si  atom  has  a  propensity  to  form 
tetrahedral  bonds.  Calculations  based  on  this  potential  reproduce 
well  the  properties  of  the  bulk  liquid  or  solid  silicon^^  and  give 
reasonable  results^®'®®  for  the  structure  of  the  (100)- 2X1  surface.  It 
fails  however  to  give  the  reconstruction  of  the  (111)  surface®®*®^  and 
has  a  mixed  success®^  in  reproducing  the  structure  and  the  energy 


of  small  clusters  given  by  ab  initio  calculations.^®  Several 
modifications  have  been  made  to  improve  the  success  of  the  potential 
for  specific  applications, 39-41  but  here  we  use  the  original  potential. 
Some  authors  believe  that  the  SW  potential  tends  to  fail  when  dealing 
with  over  or  imder  coordinated  atoms,®"^  but  that  it  is  adequate  for 
"qualitative"  studies  of  the  properties  of  the  (100)  surface. 

There  are  many  other  methods  for  calculating  the  interatomic 
forces.  A  large  variety  of  functional  forms  were  proposed  by 
Keating, Pearson  et  Takai  et.  Tersoff,^® 

Biswas  and  Hamann,^'*  Drison,^®  Brenner  and  Garrison,^® 
Chelikowski  et.  Khor  and  Das  Sarma,^®  and  Bolding  and 

Anderson.®*^  To  these  we  add  electronic  structure  methods^®  (which 
are  often  of  uncertain  accuracy)  including  those  used  simultaneously 
with  the  molecular  dynamics  calculations,  such  as  the  imbedded 
atom,®®  the  Car-Parinello®^  and  tight  binding  methods. 

Very  few  calculations  of  the  quantities  of  interest  here  have  been 
performed  with  other  methods.  The  binding  energies  for  a  single 
adsorbed  Si  atom  and  the  energy  barriers  between  the  adsorption 
sites  have  been  calculated  with  the  density  functional  method®®*®  l  and 

a  modified  Tersoff  potential.^^^'^^  Their  results  are  in  qualitative 
agreement  with  those  found  here  but  the  magnitudes  of  the  barriers 
between  sites  and  the  precise  position  of  the  binding  sites  differ.  At 
this  time  there  is  no  reason  for  preferring  one  method  over  the 
others. 

The  existing  experiments  cannot  be  used  as  a  decisive  test  of  the 
quality  of  a  given  potential,  since  there  are  as  yet  no  "zero  coverage" 
measurements  that  would  allow  determination  of  the  energy  barriers 
with  a  high  degree  of  confidence.  The  efforts  to  infer  a  diffusion 
constant  out  of  high  coverage  measurements  use  rather  sophisticated 
modeling  whose  validity  is  uncertain.  Moreover,  the  migration  at 
higher  coverages  is  rather  complicated  and  it  may  have  little 
connection  to  the  mobility  of  isolated  atoms. 

The  main  purpose  of  the  calculations  presented  here  is  to 
suggest  a  kinetic  mechanism  through  which  the  growth  patterns 
observed  on  the  Si(lOO)  surface  are  generated,  rather  than  to  evaluate 
precisely  the  binding  energies  of  all  the  adsorption  sites  and  the 


magnitude  of  the  barriers  between  them.  In  such  a  pursuit  it  is  often 
sufficient  to  find  out  which  jumps  are  the  most  frequent,  without 
having  to  know  the  precise  numerical  values  of  the  jumping  rates. 
The  fact  that  our  calculations  lead  to  growth  patterns  that  are  very 
similar  to  those  observed  experimentally  is  pleasing;  while  this 
might  indicate  that  the  SW  potential  is  chemically  reasonable,  it  does 
not  imply  that  it  is  accurate.  Further  testing  of  the  suggestions  made 
here  by  detailed  STM  measurements  is  highly  desirable. 

An  extremely  important  practical  feature  of  molecular 
dynamics  and  Monte  Carlo  calculations  is  the  need  to  make  an 
enormous  number  of  evaluations  of  the  potential  energy.  Even 
though  we  have  optimized  these  calculations  and  the  SW  potential  is 
relatively  simple,  we  are  frequently  straining  the  available  computer 
power.  A  potential  in  which  increased  accuracy  is  obtained  at  the 
expense  of  the  evaluation  time,  is  not  a  practical  option. 


a  ADSORPTION 

A  study  of  Si  adsorption  on  Si(100)-2X1  surface  should  answer 
several  questions:  (1)  What  are  the  binding  sites  for  a  newly  adsorbed 
atom?  (2)  Does  the  adsorption  of  a  single  atom  change  the  structure 
of  the  surface  layer?  (3)  How  are  the  incident  and  the  binding  energy 
of  the  incoming  atom  dissipated  when  sticking  takes  place?  (4)  Are 
there  any  low  coverage  collective  effects? 

In  many  cases,  the  deposition  of  an  atom  on  a  reconstructed 
surface  will  induce  interesting  displacements  in  the  surface  layer. 
As  more  and  more  atoms  are  deposited,  the  surface  layer  becomes  a 
sub-surface  layer.  At  some  point  during  deposition  the  former 
surface  atoms  must  unreconstruct  and  take  the  bulk  equilibrium 
positions. 

This  transition  can  take  place  in  two  different  ways.  If  may  be  a 
phase  transition  in  which  a  large  number  of  newly  deposited  atoms 
are  needed  to  induce  the  collective  rearrangement  of  a  large  group  of 
surface  atoms.  The  alternative  is  that  the  unreconstruction  of  the 
surface  is  a  local  effect  determined  mainly  by  the  interaction  between 
the  newly  arrived  atom  and  the  nearby  surface  atoms.  In  this  case 


the  unreconstruction  takes  place  on  a  first  come  first  served  basis, 
and  should  be  observable  in  the  early  stages  of  the  deposition  studied 
here. 

3.1.  Single  Atom  Adsorption  Dynamics 

In  studying  the  adsorption  dynamics  we  model  the  surface  by  a 
cluster  of  99  atoms,  distributed  in  five  layers.^^  When  molecular 
d3mamics  calculations  are  performed  14  atoms  are  propagated  by 
solving  Newton's  equations,  and  29  obey  a  Langevin  equation  using 
the  correct  potentials,  a  local  memory-less  fHction  and  a  Gaussian 
random  force  related  to  the  friction  through  the  customary 
fluctuation-dissipation  theorem.  The  rest  of  the  atoms  in  the  cluster 
are  held  fixed  and  provide  a  template  for  the  atoms  that  are  moving. 

Prior  to  the  deposition  of  a  new  atom,  the  substrate  atoms  that 
will  be  allowed  to  move  during  the  MD  simulation  are  relaxed  to  their 
equilibrium  positions  by  using  the  Metropolis  Monte  Carlo 
procedure.®®  The  relaxation  is  carried  out  at  the  temperature  of  500 
K,  which  is  close  to  570  K,  the  temperature  at  which  good  quality 
epitaxial  growth  has  been  observed.^ 

We  choose  the  center  of  a  surface  dimer  to  be  the  origin  of  the  xy- 
plane.  The  plane  of  the  centers  of  the  outermost  atoms  before  surface 
reconstruction  defines  z=0.  The  incident  Si  atom  starts  at  a  position 
(0,  0,  z)  well  above  the  surface,  where  the  interaction  energy  between 
the  incident  atom  A  and  the  surface  is  practically  zero.  The  initial 
velocity  of  A  is  perpendicular  to  the  surface  and  its  magnitude  is 
sampled  from  a  Maxwell  -  Boltzmann  distribution  at  a  specified 
temperature.  For  each  trajectory,  the  active  substrate  atoms  start 
from  a  spatial  configuration  produced  by  Monte  Carlo  sampling,  with 
velocities  sampled  from  a  Maxwell  -  Boltzmann  distribution. 

We  denote  by  T1  and  T2  in  Figs.  2(a)  and  2(c)  the  sites  occupied  by 
the  atoms  of  a  newly  deposited,  defect  free  surface  layer.  Our 
calculations  show  that  an  incident  atom  aimed  at  the  site  T1  breaks 
the  dimer  bond,  inserts  itself  into  the  dimer  and  gets  trapped  there. 
In  the  final  configuration  which  we  call  a  trimer  (see  Fig.  2(e)),  the 


two  atoms  formerly  engaged  in  the  dimer  move  to  their  equilibriiim 
bulk  positions,  with  the  adatom  residing  at  the  apex. 

This  direct  insertion  dimer  opening  mechanism  takes  place 
frequently,  and  is  insensitive  to  the  simulation  parameters  used.  We 
have  independently  varied  the  temperature  of  the  incident  atom  and 
the  surface  temperature  from  0  to  500  K,  and  the  impact  parameter 
was  sampled  to  cover  uniformly  the  square  1x1,  I  y  I  <  0.77  A  (Fig. 
2(a)).  In  all  cases  at  least  50%  of  the  trajectories  form  a  trimer  by 
direct  insertion  into  the  dimer.  The  direct  insertion  mechanism  for 
adsorption  and  dimer  opening  was  also  observed  recently  by 
Srivastava  and  Garrison  using  the  Tersoff  potential,^^  while  in 
earlier  simulations  its  importance  has  not  been  emphasized.56,57 

The  Stillinger  -  Weber  potential  gives  a  dimer  bond  strength  of 
nearly  2  eV.  This  bond  is  broken  simultaneously  with  the  formation 
of  two  new  bonds  between  the  incident  atom  and  the  atoms  that  used 
to  form  the  dimer.  Thus,  the  insertion  of  the  incident  atom  is 
exothermic. 

Initially  most  of  the  energy  released  by  the  formation  of  new 
bonds  is  localized  on  the  trimer  as  vibrational  energy,  and  it  is 
transferred  to  the  lattice  on  a  time  scale  longer  than  a  picosecond. 
After  insertion,  the  incident  atom  oscillates,  mostly  in  the  direction 
perpendicular  to  the  surface;  the  surface  atoms  initially  bound  in  the 
dimer  oscillate,  after  the  trimer  formation,  mostly  parallel  to  the 
surface. 

An  incident  atom  that  does  not  break  the  diiuer  bond  is  most 
likely  to  end  at  an  F®  site  (see  Fig.  3(a)).  It  can  get  there  in  two  ways. 
The  incident  atom  knocks  off  one  surface  atom  from  a  dimer  and 
places  it  at  the  site  F^,  as  shown  in  Fig.  3(a);  the  incident  atoms  pairs 
up  with  the  surface  atom  left  behind  (see  Fig.  3(a))  to  reform  the 
dimer.  Or,  the  incoming  atom  is  deflected  by  the  dimer  and  settles  at 
the  F®  site.  The  atom  at  the  site  F®  is  0.5  A  below  the  surface  layer, 
and  is  unable  to  open  the  neighboring  dimers  to  form  a  trimer. 
Instead,  the  two  neighboring  dimers  are  slightly  pushed  away  from 
their  2X1-(100)  equilibrium  positions. 

The  incident  atoms  aimed  at  the  site  T2  (Fig.  2(c))  have  a 


tendency  to  stick  at  the  site  F*  (  see  Fig.  3(b)).  An  atom  in  that 
position  is  unable  to  break  the  bonds  of  the  neighboring  dimers  and 
induce  the  unreconstruction  shown  in  Fig.  2(d).  It  only  manages  to 
pull  the  neighboring  dimers  closer  to  the  F*  site  by  a  distance  of  about 
0.4  A.  F*  differs  from  F*^  only  through  its  z-coordinate;  ^ 

and  the  F*  site  is,  imlike  F®,  above  the  surface  but  below  the  apex  of  a 
trimer. 


(a)  (b)  (c)  (d)  (e) 


Fig.  2.  T1  and  T2  indicate  the  location  of  the  deposited  atoms  after  the 
deposition  of  a  monolayer,  (b)  and  (d)  show  possible  surface  atom 
displacements  caused  by  the  adsorption.  Only  (b)  takes  place  during 
the  formation  of  the  trimer  shown  in  (e).  The  solid  circles  represent 
adatoms;  the  stripped  are  the  substrate's  surface  atoms. 


Fig.  3.  Adsorption  processes  leading  to  the  occupation  of  the  bridge 
sites  between  two  dimer  rows.  F®  in  (a)  is  below  the  surface  layer;  F* 
in  (b)  is  above  the  surface  layer. 


3^2.  Collective  Effects 


As  noted  above,  an  atom  adsorbed  half  way  between  two  dimers 
is  either  below  the  surface  (if  adsorbed  at  F®)  or  above  the  surface  (at 

F*)  but  below  the  new  layer.  In  addition,  such  an  adatom  is  unable  to 
pull  the  nearby  surface  atoms  out  of  position  and  form  a  trimer  with 
them.  Since  trimerization  must  take  place  by  the  time  the  deposition 
of  a  monolayer  is  completed,  we  infer  that  the  atoms  that  happen  to 
stick  to  the  F*  or  the  F®  positions  can  only  induce  trimerization  if 
assisted  by  atoms  deposited  subsequently. 

In  searching  for  such  a  collective  effect,  we  have  placed  an  atom 
A1  at  F®  and  performed  a  large  number  of  MC  moves  to  make  sure 
that  all  surface  atoms  were  in  positions  compatible  with  thermal 
equilibrium.  This  means  that  they  were  located  at  positions 
generated  by  importance  sampling,  not  at  the  mean  equilibrixim 
positions. 

After  the  equilibration  is  completed  we  aim  an  atom  A2  at  T1 
with  a  thermal  velocity  perpendicular  to  the  surface,  as  shown  in 
Fig.  4(a).  Many  of  these  trajectories  lead  to  the  trapping  of  A2  at  the 
center  of  the  dimer,  labeled  as  the  B  site.  This  breaks  the  dimer  bond 
S3  -  S4,  and  the  trimers  S2-  A1  -  S3  and  S3  -  A2  -  S4  (see  Fig.  4(b))  are 
formed.  The  atoms  in  these  trimers  relax  to  their  bxdk-like  positions; 
compared  to  the  initial  configuration,  the  distance  S3  -  S4  is  longer, 
S2  •  S3  is  shorter,  and  A1  is  lifted  up.  Notice  that  the  dimer  SI  >  S2  is 
unaffected:  Si  follows  S2  and  maintains  its  dimer  bond  with  it. 

Similar  collective  effects 'are  observed  if  we  initially  place  an 
atom  in  equilibrium  at  F*  and  send  a  second  atom  at  B. 

These  collective  effects  resemble  those  observed  by  Srivastava, 

Garrison,  and  Brenner®^  in  simulations  using  the  Tersoff  potential. 
However,  in  their  simulations  the  adatoms  deposited  earlier  first 
saturate  the  two  dangling  bonds  of  a  dimer  (the  A  sites  in  our 
notation  in  Sec.  3.3),  then  a  third  adatom  breaks  the  dimer  bond  and 

inserts  in  it.®^  Our  results  suggest  that  such  a  mechanism  is 
improbable  at  the  low  deposition  rates  used  in  epitaxy,  since  an 
adatom  deposited  in  an  A  site  will  hop  to  an  F  site  before  the  arrival  of 

another  adatom.^® 


Fig.  4.  Collective  effect:  The  arrival  of  A2  helps  A1  to  take  the 
position  it  will  have  when  the  layer  is  completed. 


3.3.  Ibe  Distribution  of  the  Adsorption  Sites 

Next  we  deposit  the  Si  atoms  uniformly  onto  the  iinit  cell  of  the 
reconstructed  Si(100)-(2X1)  surface.  Our  purpose  is  to  determine  all 
the  adsorption  sites,  and  the  site  occupation  probabilities  in  the  early 

stages  of  adsorption.^^ 

At  zero  temperatures  the  adsorption  sites  are  local  potential 
energy  minima.  We  deposit  a  particle  on  the  surface  and  use  an 
energy  minimization  program  that  moves  the  surface  atoms  and  the 
adsorbed  one  until  all  the  local  minima  are  found.  The  same 
procedure  is  used  to  move  the  adsorbed  atom  from  one  adsorption  site 
to  another  along  the  minimxim  energy  path;  this  gives  the  energy 
barriers  between  the  sites.  The  results  are  shown  in  Fig.  5  and  Table 


I.  There  are  thirteen  local  minima  within  the  surface  unit  cell,  of 

which  six,  labelled  A,  B,  C,  D,  E,  and  F,  are  distinct.®®  The  others  are 
eqmvalent,  by  symmetry,  to  one  of  the  six.  The  reference  energy  for 
the  binding  energies  at  the  local  minima  and  the  energy  barriers 
separating  them  is  that  of  the  (2X1)  surface  with  the  adsorbate  at 
infinity. 

Now  we  examine  the  adsorption  of  the  Si  atoms  on  the  surface 
sites  described  above.  In  these  calculations  an  incident  Si  atom  starts 
at  the  point  (x,  y,  3.5  A)  where  the  atom  >  surface  interaction  is  zero. 
The  X"  and  y-coordinates  cover  uniformly  the  unit  cell  defined  by  the 
rectangle  ixl<1.0  a,  lyl^.O  a,  where  a=3.84  A  is  the  surface  lattice 
constant  before  reconstruction  (Pig.  5). 
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Fig.  5.  The  adsorption  sites 
within  the  surface  unit  cell. 


Fig.  6.  A  snapshot  of  1000 
adatoms  at  t  -  0.9  ps. 


labkl 


Tlie  coordinates  (x,y,z)  of  the  binding  sites,  the  binding 
energy  Ej,,  and  the  barrier  height  Vj,  between  two  neighboring  sites,  s 
is  the  total  area  of  a  given  type  of  adsorption  site  within  the  dotted 
rectangle  in  Fig.  5. 


site 

X  (A) 

y(A) 

z  (A) 

Eb(eV) 

Vb  (eV) 

8(A2) 

A 

0.0 

2.06 

2.037 

2.162 

A->B 

0.159 

9.59 

A^F 

0.323 

A-^D 

0.200 

B 

0.0 

0.0 

1.359 

2.696 

B-^A 

0.693 

3.69 

B-^E 

0.288 

C 

1.92 

0.0 

1.508 

2.639 

C-»D 

0.176 

1.10 

C  E 

0.159 

D 

1.92 

1.73 

1.704 

2.635 

D-^C 

0.190 

8.66 

A 

0.692 

D-»F 

0.759 

E 

1.25 

0.0 

1.076 

2.613 

E->B 

0205 

1.02 

E-»C 

0.133 

F 

0.0 

3.84 

1.128 

3.140 

F->  A 

1.301 

5.90 

F-»  F* 

0.898 

F->D 

1245 

As  one  may  expect,  the  atoms  are  drawn  towards  a  few  sites 
within  the  unit  cell,  even  though  the  initial  trajectories  cover 
uniformly  the  whole  area.  What  is  perhaps  unexpected  is  the 
population  inversion  observed  in  the  earliest  stages  of  adsorption:  the 
incident  atom  is  adsorbed  with  high  probability  at  the  sites  A,  B,  C 
and  D  in  Fig.  5,  but  not  at  the  F  site,  even  though  the  latter  is  the  site 


with  the  largest  binding  energy;  the  site  A  having  the  smallest 
binding  energy  attracts  the  largest  fraction  of  the  incident  atoms. 

In  Fig.  6,  we  show  a  snapshot  of  the  x  and  y  coordinates  of  1000 
particles  at  0.9  ps.  It  is  clear  that  at  this  early  time  the  trajectories 
have  already  clustered  up  around  the  sites  A,  B,  C  and  D.  By 
following  the  particles  for  a  while  longer  we  observe  (see  Fig  7)  a  slow 
depopulation  of  A  accompanied  by  a  population  gain  in  B,  which  goes 
on  for  five  picoseconds.  This  means  that  on  this  time  scale  the  kinetic 
energy  of  some  of  the  particles  in  A  still  exceeds  the  small  energy 
barrier  (  0.159  eV)  between  A  and  B. 

In  Fig.  7  we  show  the  time  evolution  of  the  fraction  of  trajectories 
whose  X  and  y  coordinates  fall  within  the  squares  that  define  the 
adsorption  sites  A,  B,  C,  and  D.  The  initial  population  is  proportional 
to  the  area  associated  with  a  given  site  (see  the  last  column  of  Table 
I).  After  about  0.2  picoseconds  the  atom  -  surface  interactions  steer 
th<»  incident  particle  towards  the  site  A,  away  from  F,  D  and  E. 
Practically  no  particle  manages  to  adsorb  at  F.  The  site  occupation  is 
pretty  much  established  in  about  0.6  ps,  and  the  site  A  captures  the 
largest  fraction  of  particles. 
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Fig.  7.  The  occupation  of  the  surface  sites  at  the  early  times  of 
adsorption  for  300  adatoms. 


4.  THE  MIGRATION  OF  AN  ISOLATED  ADATOM  ON  THE 
SURFACE 


4.L  Hie  Site>to>SiteHoppii]^  Rates 


At  the  temperatures  of  interest  here  the  migration  of  an 
adsorbed  atom  on  the  surface  takes  place  through  thermally  activated 
site  to  site  jumps.  The  rates  of  such  jumps  can  be  crudely  estimated 
by  using  the  Arrhenius  formiila  with  a  pre-exponential  of  10^^  sec'^ 
and  the  energy  barriers  given  in  Table  I.  Methods  for  exact 
calculations  of  the  jumping  rates  are  available®®*®^  but  will  not  be 
used  here  since  our  purpose  is  to  gain  a  qualitative  understanding  of 
growth  phenomena.  In  Fig.  8  we  show  the  site  to  site  hopping  rates 
at  T  =  500  K,  divided  by  the  rate  of  the  E  ->  C  jump. 


Fig.  8.  The  site  to  site  hopping  rates  at  500  K. 


The  motion  in  the  direction  perpendicular  to  the  dimer  rows  (i.e. 
B  A,  D  F,  F  ->  A  or  F  D)  is  extremely  slow.  Furthermore,  an 
atom  adsorbed  at  B  can  move  easily  along  the  line  B  E  C  E* 
B*  with  likely  side  trips  from  C  to  D  and  back.  Since  there  is 
essentially  no  fast  pathway  for  an  atom  to  exit  the  chain  BEC(D)E*B*, 
the  early  motion  after  adsorption  is  practically  one-dimensional.  On 
a  longer  time  scale  the  atom  has  a  fair  chance  to  jump  into  the  site  F 
where  it  will  rest  for  a  long  time.  On  this  time  scale  the  sites  F  act  as 
traps.  At  an  even  longer  time  scale  the  atom  can  move  to  F*,  or,  on 
the  longest  time  scale,  come  back  to  D  or  A. 

We  note  that  the  anisotropy  of  diffusion  on  Si(lOO)  has  been 
suggested  by  Legally  et  al.  to  explain  their  STM  results.5.l3  Lu  and 
Metiu®^^  have  argued  that  such  anisotropy  should  be  expected  for  all 
(100)  faces  of  semiconductors  having  diamond  or  zincblende 
structures,  since  the  barrier  along  the  dimer  rows  is  lowered  by  the 
fact  that  the  breaking  of  the  bonds  to  the  initial  site  is  concerted  with 
the  formation  of  bonds  with  the  final  site;  no  such  lowering  is  present 
when  the  atom  is  moved  in  a  direction  perpendicular  to  the  dimer 
rows  and  the  energy  barriers  for  that  kind  of  displacement  ought  to  be 
higher. 

4.2.  Diffusion  on  a  Lattice  with  Many  Sites 

Most  simulations  of  epitaxial  growth  have  used  models  in  which 
a  particle  jumps  randomly  on  a  square  lattice.  It  is  well  known  that 
for  such  a  system  the  probability  of  finding  the  particle  at  a  given  site 
is  described  by  a  diffusion  equation  and  the  diffusion  coefficient  is 
proportional  to  the  jumping  rate.  The  temperature  dependence  of  the 
diffusion  coefficient  is  the  same  as  that  of  the  jumping  rate  and  the 
latter  is  given  by  the  Arrhenius  formula.  Our  calculations  suggest 
that  reality  is  more  complicated.  On  a  short  time  scale  the  random 
motion  takes  place  along  a  string  with  several  sites  and  the  particle 
can  make  side  trips  to  the  site  D  and  back  with  a  fairly  high 
frequency. 

If  such  motion  is  diffusional,  the  diffusion  coefficiei  .c  must 


depend  on  all  the  site  to  site  rate  constants.  Since  the  Arrhenius 
formula  applies  to  the  jumping  rate,  not  to  the  diffusion  coefficient, 
the  relationship  between  the  diffusion  coefficient  and  the  temperature 
can  be  rather  complicated.  If  the  relationship  happens  to  be  of  an 
Arrhenius  form,  the  activation  energy  and  the  pre-exponential  will 
not  have  a  simple  meaning.  On  the  time  scale  on  which  the  atom  is 
likely  to  be  trapped  into  an  F  site  the  motion  is  not  diffusional;  i.e.  the 
mean  square  displacement  is  not  proportional  to  time. 

To  understand  the  nature  of  the  adatom  migration  we  have 
performed  stochastic  kinetic  simulations  to  calculate  the  mean 
square  displacements  for  three  models:  One  in  which  both  the  side 
trips  to  D  and  the  access  to  the  traps  are  forbidden;  one  in  which  only 
the  traps  are  inaccessible;  and  one  in  which  all  processes  are 
allowed.^® 

The  stochastic  kinetic  simulation  moves  a  Si  atom  along  the 
sites  shown  in  Fig.  8  with  probabilities  proportional  to  the  ratio 
between  the  rate  constant  of  the  appropriate  jump  and  a  reference 
rate  k^  We  use  for  k^  the  rate  of  the  jump  E-»C  at  600  °C,  which  is  the 
fastest.  The  E->C  jump  is  performed  whenever  attempted.  The  time 
(ky  )*^  provides  the  time  scale  for  the  simulation.  The  time  in  the 
plots  in  this  subsection  is  k^t,  where  t  is  the  clock  time.  The  time  unit 
(kj.  )'^  happens  to  be  close  to  a  picosecond.  The  averages  are 
performed  over  10^  independent  nms. 

4J2.1.  The  Random  Walk  on  the  Fast  String 

To  understand  the  migration  process  we  examine  first  the 
random  walk  on  the  fast  string  ..BECE*B*..  when  the  jumps  to  A,  D 
and  F  are  forbidden.  The  highest  barrier  along  this  string  is  0.29  eV 
(for  B->E),  and  the  lowest  is  0.13  eV  (for  E->C).  The  time  dependence 
of  the  mean  square  displacement  (x^)  obtained  in  the  stochastic 
kinetic  simulations  at  400,  600,  and  800  ®C  are  shown  in  Fig.  9.  The 
plots  are  linear  in  t,  indicating  that  the  motion  is  diffusional;  the 
proportionality  constant  is  twice  the  diffusion  coefficient.^® 
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Fig.  9.  Mean  square  displacements  as  a  function  of  time  at  400,  600 
and  800  ®C,  for  the  three  models. 


The  temperature  dependence  of  the  dif^sion  constants  is  plotted 
in  Fig.  10,  for  a  temperature  range  from  200  to  800  °C.  This  curve  has 
an  Arrhenius  form  with  an  activation  energy  of  0.24  eV  and  a  pre¬ 
exponential  of  2.36  /x.  Here  x=(ky)'^  and  the  reference  rate  = 

1.71x10^^  sec'^.  We  can  define  an  effective  j\imping  rate  kg  by  using 
Dg=  a^  kg  .where  a  is  the  B  to  B*  distance.  The  pre-exponential  Ag  in 
the  effective  jumping  rate  kg  is  Ag  =  2.36  k^/a^  =  2.76x10^^  sec*^.  We 
see  here  a  "compensation  effect":  while  all  the  rates  in  the  model 
have  pre-exponentials  of  10^^  sec'^  the  pre-exponential  of  the  effective 
rate  is  almost  two  orders  of  magnitude  smaller.  The  activation 
energy  of  the  diffusion  coefficient  is  also  smaller  than  that  of  the 
highest  barrier  along  the  fast  string. 


Fig.  10.  Diffusion  coefficient  for  two  different  models  of  multisite  one¬ 
dimensional  diffusion:  without  and  with  side  trips,  respectively. 


To  understand  these  restilts  we  have  examined  a  simpler  model 
in  which  an  atom  jumps  along  a  string  which  has  only  two  sites  per 
tmit  cell  (Fig.  11).  The  rate  equation  for  this  system  is 

=  k2iP2(a-l;t)  +  k2iP2(a;t)  -  2  ki2Pi(a;t)  (4.1a) 

dt 

=  ki2Pi(a+l;t)  +  ki2Pi(a;t)  -  2  k2iP2(a;t) .  (4.1b) 

dt 

Pj(a;t)  is  the  probability  that  the  atom  is  at  a  site  of  type  i  (i=:l  or  2)  in 
cell  a,  at  time  t,  and  k^j  is  the  rate  of  the  i-»j  jump. 


Fig.  11,  A  schematic  representation  of  an  one  dimensional  difhision 
channel  with  two  inequvalent  sites. 


A  simple  analysis^®  shows  that  the  probability  P(a;t)= 
Pl(oc;t)+P2(a;t)  that  the  particle  is  in  the  cell  a  at  time  t,  satisfies  a 
diffusion  equation  with  the  diffusion  constant 


D=  1  { }  a2. 
2  ki2+k2i 


(4.2) 


Clearly  the  diffusion  coefficient  does  not  have  a  simple 
Arrhenius  temperature  dependence  if  each  of  the  two  rate  constants 


k]^2  ^21  Nevertheless,  as  seen  in  the  simulations  reported 

above  the  temperature  dependence  is  often  fitted  well  by  the 
Arrhenius  formula.  To  find  the  meaning  of  the  pre-exponential  and 
the  activation  energy  in  such  fits  we  expand  In  D  in  powers  of  l/(k3T): 


2  fi2+f21  fl2  +  f21  kfiT 


r  fl2  f21  (Vl2  -  V2l)^  i  /_ 

2(fl2+f2l)^  ^bT 


f  +  0( 


bT 


(4.3) 


Here  and  Vy  are  the  pre-exponentials  and  the  activation  energies 
in  the  rate  constants  k^j.  The  pre-exponential  in  the  diffusion 
coefficient  is  a^  ^12^21  ^  activation  energy  is  (fi2V2i 

■*■^21  ^12  )/(fi2‘*'^2i^-  Arrhenius  formula  is  valid  when  the  third 
term  on  the  right  hand  side  is  negligible  compared  to  the  first  two. 
Similar  relationships  exist  for  three  and  four  site  systems.^® 

4JL2.  Diffusion  Along  a  String  with  a  Side  Trip 


Qiialitatively,  it  is  not  difficiilt  to  see  that  the  side  trips  from  the 
C  to  the  D  sites  and  back  will  delay  the  progress  of  the  atom  along  the 
fast  string,  thus  leading  to  a  smaller  diffusion  coefficient.  To  find 
how  important  this  effect  is  we  performed  stochastic  kinetic 
simulations  with  a  model  in  which  these  side  trips  are  allowed,  but 
the  jumps  to  A  and  F  are  still  forbidden.  The  time  dependence  of  the 
mean  square  displacement  along  the  fast  string  is  shown  in  Fig.  9  for 
400, 600,  and  800  ®C.  The  proportionality  to  t  means  that  the  motion  is 
diffiisional.  The  temperature  dependence  of  the  diffusion  coefficient 
obtained  from  the  mean  square  displacement  satisfies  the  Arrhenius 
equation  (see  Fig.  10)  and  a  fit  of  the  results  leads  to 

DCD  =  1.23  exp{  -  ^ }  (A2  k^ )  (4.4) 

Kb  1 


Allowing  the  side  trips  did  not  change  the  effective  barrier,  but 


diminished  the  pre-exponential  factor  from  2.36  to  1.23  A^k,.. 

The  analysis  of  a  simplified  model  of  random  jumps  along  a 
chain  having  two  sites  in  the  unit  cell  along  the  chain  eind  one  site  for 
a  side  trip  (see  Fig.  12)  gives  the  diffusion  coefficient 

D  _ ]  (4  5) 

2  (ki2+k2l)  +  ki2k23/k32 

The  definition  of  the  rate  constants  k^j  is  given  in  Fig.  12.  Note  that 
this  equation  is  not  valid  if  k32  becomes  zero  because  in  that  case  the 
site  3  acts  as  a  trap  and  the  motion  is  no  longer  diffusional.  In  the 
other  extreme  where  k32  is  very  large  or  k23  is  very  small,  this 
equation  reproduces  the  diffusion  coefficient  for  the  chain  without  a 
side  trip.  A  chain  with  a  side  trip  is  identical  to  a  chain  without  a 
side  trip  in  which  the  lifetime  at  the  site  from  which  the  side  trip 
originates  is  increased. 


a-1  a  a-(-l 


Fig.  12.  A  schematic  representation  of  a  model  in  which  the  particle 
can  jump  along  an  one  dimensional  lattice  having  two  sites  (the 
square  and  the  circle)  and  it  can  also  take  a  side  trip  from  the  site  2  to 
the  site  3  ( the  ellipse). 


42.3.  Migration  Along  the  Fast  String  with  Side  Trips  and  Traps 

Now  we  present  the  resvilts  of  sirntdations  in  which  both  the  side 
trips  and  the  traps  are  included.  The  time  dependence  of  the  mean 
square  displacements  <x^>  is  shown  in  Fig.  9.  The  curves  for  <x^> 
rise  linearly  at  early  times,  then  bend  and  flatten  out  at  longer  times. 
This  happens  because  there  is  a  small  probability  that  a  Si  atom 
visiting  the  site  D  will  jump  to  the  site  F  and  remain  trapped  there. 
When  this  happens  the  atom  no  longer  contributes  to  the  growth  of 
<x^>.  The  time  dependence  of  the  mean  square  displacement  varies 
smoothly  from  linear  to  constant  because  different  atoms  are  trapped 
at  different  times.  At  a  time  scale  longer  than  the  lifetime  at  the  site 
F  the  visiting  atom  leaves  the  trap  by  a  hop  parallel  to  the  dimer  rows 
into  another  trap  F*;  the  barrier  for  this  move  is  0.9  eV.  The 
probability  to  hop  back  onto  the  fast  string  is  rather  small,  with  a 
barrier  of  1.25  eV.  The  migration  on  this  longer  time  scale  is  not 
investigated  here. 

Notice  the  mean  migration  distance  at  800  °C  is  shorter  than 
that  at  600  ®C.  This  happens  because  at  the  higher  temperature  the 
trapping  rate  is  larger  and  the  mean  square  displacement  starts 
deviating -from  the  linear  behavior  at  an  earlier  time. 

The  observed  behavior  can  be  understood  by  examining  diffusion 
along  an  one-dimensional  string  with  a  small  probability  that  the 
particle  leaves  the  string  and  gets  trapped.  To  describe  this  system 
we  must  use  a  random  walk  equation  with  a  rate  term  to  jump  into  a 
trap,  and  no  possibility  of  coming  back.  The  mean  square 
displacement  <x^>  for  such  a  model  is 

Here  D  is  the  diffusion  coefflcient  on  the  string  and  s  is  the  rate 
constant  to  fall  in  the  trap.  To  test  whether  this  equation  describes 
well  the  behavior  of  <x^>  observed  in  the  simulation  we  calciilate  the 
mean  square  displacement  for  D  =  0.07  k,.,  which  is  the  diffusion 

constant  of  the  one  dimensional  motion  with  side  trips,  and  s=10*^. 


which  is  estimated  from  the  rate  of  hopping  from  C  to  D  and  then  D  to 
P.  The  empty  circles  in  Fig.  13  show  the  results  given  by  Eq.  (4.6) 
while  the  solid  line  shows  the  results  of  the  simulation.  The 
agreement  is  excellent. 


Fig.  13.  Time  dependence  of  the  mean  square  displacement  for  the 
model  allowing  side  trips  to  D  and  visits  to  F. 

5.  ADATOM-ADATOM  INTERACTIONS  AND  THE  PATHWAYS 
FOR  DIMER  STRING  GROWTH 

The  STM  measurements show  that  at  low  coverage  the 
silicon  atoms  adsorbed  on  the  Si(100)-2X1  surface  form  string  like 
islands  perpendicular  to  the  dimer  rows  on  the  surface.  These  island 
shapes  are  highly  unusual.  The  reason  why  the  silicon  atoms  prefer 
these  shapes  and  the  mechanism  of  their  formation  are  not 
understood.  The  large  number  of  strings  indicates  that  many 
nucleation  centers  are  produced  in  the  early  stage  of  the  deposition. 
It  is  not  clear  what  these  nucleation  centers  are.  Our  studies 
reported  in  the  previous  section  indicate  that  the  F  sites,  located 


between  the  substrate's  dimer  rows,  are  very  slow  to  fill;  the  site  to 
site  jumps  leading  to  them  have  low  rates  compared  to  those  of  the  B 
sites.  To  form  the  strings  observed  by  STM  the  migrating  atoms  fill 
as  many  F  sites  as  B  sites.  This  is  rather  surprising  and  reqmres  an 
explanation.  In  particular  we  need  to  know  why  there  are  so  few  F 
site  vacancies  in  the  islands.  Finally,  the  work  reported  in  the 
previous  section  indicates  that  the  diffusion  of  the  single  atoms  on  the 
surface  is  highly  anisotropic.  Since  the  strings  are  perpendicular  to 
the  direction  of  fast  diffusion,  most  of  the  single  atoms  moving  on  the 
surface  will  reach  the  side  of  an  existing  string;  given  this,  it  is 
difficult  to  understand  how  the  string  is  getting  longer  and  not  wider. 
Below  we  provide  an  answer  to  these  questions. 

In  many  surface  science  systems  it  is  assumed,  with  good 
reason,  that  the  diffusion  coefficient  at  high  coverage  is  not  radically 
different  from  that  of  a  single  particle.  The  change  caused  by 
increased  coverage  is  usually  calculated  by  a  mean  field  model  in 
which  the  change  depends  on  the  mean  coverage  only.  This  is  not 
true  for  silicon,  for  which  the  interactions  are  not  only  very  strong  but 
are  also  spatially  directed.  The  rate  to  jump  from  site  to  site  depends 
not  only  on  the  number  of  neighbors  surroxmding  the  sites^'^^  but 
also  on  their  precise  location.  It  is  useful  to  think  of  the  two 
dimensional  clusters  formed  by  the  adsorbed  atoms  as  distinct 
isomers:  clusters  with  the  same  number  of  atoms  have,  in  most 
cases,  substantially  different  properties  if  their  "structure"  differs. 
The  isomerization,  dissociation  and  atom  capture  rates  differ  from 
isomer  to  isomer.  To  understand  how  the  islands  grow  and  change 
shape  we  must  consider  the  properties  of  a  very  large  number  of 
clusters.  We  discuss  here  only  those  that  provide  the  fastest  steps^S 
leading  to  the  strings  observed  experimentally.  The  rates  are 
estimated  by  calculating  the  energy  barriers  for  each  step  and 
assuming  that  all  rate  constants  have  the  same  pre-exponential. 

We  consider  configuration  changes  achieved  by  moving  one 
atom  from  a  site  i,  of  the  initial  configuration  to  a  site  j  of  the  final 
one.  Occasionally,  we  have  also  investigated  changes  achieved  by 
moving  dimers,  when  our  intuition  suggested  that  they  might  be 
rapid;  we  have  found  dimer  motion  too  slow  to  be  an  important  factor 
in  island  growth.  The  potential  energy  barrier  to  go  from  a  site  to 


another  site  Sj  is  determined  as  follows.  We  first  place  the  jumping 
atom  at  S|.  Then  the  atomic  positions  of  all  the  atoms  (the  one  whose 
motion  is  studied,  the  adsorbed  atoms  neighboring  Sj  and  Sj,  and  the 
pertinent  substrate  atoms)  are  relaxed  to  minimize  the  total  energy  of 
the  system.  Next,  we  change  the  position  of  the  atom  whose 
migration  we  want  to  study,  by  a  small  displacement  along  the 
surface  towards  the  final  site,  by  changing,  for  example,  its  x 
coordinate.  Then  we  keep  x  fixed  at  this  new  value  and  vary  the 
coordinates  of  all  the  other  atoms,  as  well  as  the  y  and  z  coordinates 
of  the  atom  of  interest,  until  we  find  the  minimum  energy.  This 
procedure  is  repeated  until  the  atom  being  displaced  reaches  the  final 
site. 

The  interaction  with  the  substrate  and  the  neighboring  atoms  is 
so  strong  that  it  is  not  possible  to  describe  the  reaction  path  by 
specifying  only  the  displacement  of  the  jumping  atom.  As  the 
jumping  atom  moves  towards  the  final  site,  the  other  atoms  can  also 
undergo  substantial  displacements.  There  are  in  fact  many 
pathways  by  which  an  atom  goes  from  Sj  to  Sj,  that  differ  through  the 
way  the  other  atoms  are  moving,  and  have  different  energy  barriers. 
The  one  with  the  lowest  barrier  defines  the  preferred  jumping  path. 

The  substrate  is  modeled  by  a  cluster  of  495  atoms,  distributed  in 
five  layers.  Of  these,  213  atoms  are  allowed  to  relax  to  give  the  lowest 
potential  energy  for  the  whole  system;  the  rest  are  fixed  and  provide  a 
template.^® 

At  this  point  it  is  useful  to  remember  that  an  atom  landing  on  a 
Si(100)-2X1  surface  moves  rapidly  along  the  dimer  rows,  and  the  rate 
to  jtunp  to  an  F  site  is  low.  Nevertheless,  at  the  low  deposition  rates 
typical  to  epitaxial  growth  experiments,  the  atoms  deposited  earliest 
are  likely  to  reach  an  F  trap  before  having  a  chance  to  encounter 
another  adsorbed  atom.  The  most  likely  path  to  jump  into  an  F  site 
from  the  fast  channel  is  via  C->D-»F,  with  barriers  of  0.15  and  0.68 
eV,  respectively.  The  jump  of  an  atom  from  an  F  site  back  on  the  fast 
string  is  very  improbable.  The  F-^F*  jump  also  has  a  large  barrier 
of  0.9  eV  (see  Fig.  8).  For  this  reason  we  call  the  set  of  sites  F,  F*,  etc. 
located  between  the  substrate  dimer  rows,  the  slow  string. 

Atoms  deposited  later  move  along  the  fast  strings  and  either  fall 


into  unoccupied  F  sites  or  bind  to  the  occupied  ones  forming  pairs, 
with  one  atom  in  F  and  the  other  in  B  (Fig.  14(b))  or  in  E  (Fig.  14(c)). 
The  pairs  are  easily  formed:  the  highest  barrier  for  the  move  from  B' 
to  B  (Fig.  14(a))  is  0.24  eV.  Once  formed,  they  are  fairly  persistent: 
the  atom  at  F  will  not  move,  and  the  barrier  to  move  from  E  to  E'  (  see 
Fig.  14(c))  is  =  0-73  eV.  Thus,  an  atom  at  F  pins  down  an  atom 

running  along  the  nearby  fast  string  localizing  it  in  either  the  B  or 
the  E  site  nearby;  the  barriers  =  0.23  eV  and  =  0.20  eV 

are  low  and  the  atom  jumps  back  and  forth  between  B  and  E. 
Therefore,  the  nucleation  centers  for  the  strings  are  the  early  atoms 
which  reach  the  F  sites.  There  are  many  such  nucleation  centers 
because  the  deposition  rate  is  very  low  and  this  gives  the  early  atoms 
a  chance  to  reach  the  F  sites,  get  stuck  there  and  trap  another  atom 
running  along  the  fast  strings. 


Fig.  14.  An  adatom  trapped  in  the  site  F  pins  down  an  atom  running 
on  the  nearby  fast  string. 


We  have  also  examined  the  mobility  of  a  dimer  located  on  the 
fast  or  the  slow  string  to  see  if  dimer  diffusion  plays  an  important 
role  in  the  growth  process.  We  find  that  once  a  dimer  is  formed  it  is 
essentially  immobile. 

As  more  atoms  are  deposited  on  the  surface  they  move  rapidly 
along  the  dimer  rows  to  form  the  triplets  shown  in  Fig.  15.  In  Fig. 
15(a),  the  three  atoms  shown  are  located  at  sites  corresponding  to  the 
growth  pattern  observed  in  STM.  The  atoms  at  and  B2  form  a 
stable  dimer.  The  presence  of  atoms  at  B^  and  B2  reduces  the  barrier 
for  the  trapped  atom  to  move  from  F^  to  F2  but  it  will  not  break  away 
from  the  trimer.  In  Fig.  15(b),  the  two  atoms  located  at  the  sites 
and  E2  are  in  the  wrong  place  (i.e.  not  where  observed 
experimentally,  in  the  final  structure). 


(a)  .  (b) 


Fig.  15.  The  formation  of  triplets:  (a)  is  stable,  (b)  is  a  precursor  state 
for  the  process  shown  in  Fig.  16. 


It  turns  out  that  the  wrong  configuration  shown  in  Fig.  150))  is 
an  important  precursor  state  for  the  occupation  of  F2  and  the 
formation  of  a  dimer  in  the  slow  string.  The  process  is  shown  in  Fig. 
16.  If  we  pull  the  atom  2  in  Fig.  16(a)  as  indicated  by  the  arrow,  the 
atom  1  will  follow  for  a  while;  then,  the  1-2  bond  breaks  and  2  ends  at 
an  P  site  (Fg  in  Fig.  16(b))  and  1  in  an  E  site  (E2  in  Fig.  16(b)).  The 
barrier  for  this  move  is  0,59  eV.  This  barrier  is  lower  than  the 
maximum  barrier  of  0.76  eV  for  the  single  atom  C  D  -»  F  hopping, 
partly  because  as  the  atom  2  breaks  the  1-2  bond  and  leaves  the  atom 
1  behind,  it  forms  a  new  bond  with  the  atom  at  F^^. 

Since  in  the  low  coverage  limit  the  rate  to  reach  the  F  sites  is 
very  small,  many  such  sites  will  be  left  vacant  unless  the  presence  of 
other  particles  help  fill  the  F  sites.  The  move  shown  in  Fig.  16  is 
important  not  only  because  it  repairs  an  error  (i.e  allows  the  pair 
located  at  the  wrong  sites  and  E2  to  move  to  correct  sites)  but  also 
because  it  provides  a  mechanism  for  populating  a  site  F  and  form  a 
stable  dimer  there. 


If  a  fourth  atom  is  brought  close  to  the  triplet  configuration 
shown  in  Fig.  15(a),  we  find  two  other  moves  that  can  efficiently 
populate  the  F  site  (Fig,  17).  In  Fig.  17(a)  the  atom  at  E  can  easily 
jump  into  the  site  F2  and  form  a  dimer  with  the  atom  at  Fi  (Fig. 
17(b)).  The  highest  barrier  along  this  path  is  0.44  eV.  An  adatom  at 
C  (Fig.  17(c))  can  also  hop  into  F2  and  form  a  dimer  with  F^  (Fig. 
17(b));  the  barrier  for  this  hop  is  0.57  eV. 

The  jumps  discussed  above  provide  a  kinetic  pathways  for  the 
growth  of  a  dimer  string  of  length  two  (  Fig.  17(b) ).  All  other  jumps 
that  we  have  tried  are  too  slow  to  alter  this  pathway  significantly. 
This  very  short  string  can  be  considered  a  nucleus  for  the  growth  of  a 
longer  string.  We  investigate  next  how  this  growth  takes  place.  To 
begin  with,  this  small  cluster  will  trap  atoms  moving  on  the  fast 
string,  at  the  sites  B3  and  B4.  Further  growth  of  the  string  takes 
place  by  filling  the  F  sites  at  the  end  of  this  "nucleus"  by  one  of  the 
mechanisms  discussed  above. 
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Pig.  17.  (a)  and  (c)  show  two  more  mechanisms  by  which  the  F  sites 
are  populated.  Both  moves  lead  to  the  configuration  (b). 


We  must  now  address  what  is  perhaps  the  most  puzzling 
feature  of  the  observed  growth  patterns.  The  string  grows  in  the 
direction  perpendicular  to  the  dimer  rows,  even  though  the  free 
atoms  moving  along  the  surface  move  along  the  dimer  rows  and 
reach  the  string  on  its  side.  This  should  lead  to  a  string  widening 
rather  than  elongation.  Therefore  we  must  find  a  mechanism  for  the 
transport  of  atoms  from  the  side  to  the  end  of  the  string. 

We  have  examined  two  possibilities.  Mo  et  suggested  that 
perhaps  the  probability  that  an  atom  sticks  at  the  side  of  the  string  is 
much  smaller  than  that  for  sticking  at  the  end.  Our  energy 
calculations  show  that  atoms  reaching  the  end  of  a  string  are  tightly 
botmd  there.  The  binding  is  much  stronger  if  the  string  ends  with  an 
F  site  than  with  a  B  site.  This  is  in  agreement  with  the  suggestion  by 
Mo  et  But,  an  atom  can  easily  reach  the  side  of  the  string  (e.g. 
the  atom  1  in  Fig.  18(a))  and  is  strongly  bound  there.  The  energy 
barrier  to  jump  away  from  the  string  is  about  2  eV.  This  is  contrary 
to  the  assumption  made  by  Mo  et 

The  other  possibility  that  we  examined  is  related  to  recent  work 
by  Lu,  Petroff  and  Metiu,®  who  have  shown  that  a  square  lattice 
model  with  anisotropic  interactions  between  the  adsorbed  particles 
leads  to  the  formation  of  elongated  islands.  In  their  model  these 
interactions  facilitate  the  migration  along  the  side  of  the  string: 
particles  reaching  the  side  of  the  string  move  along  it  until  they  reach 
the  head  and  get  stuck  there.  The  present  calculations  do  not  support 
this  mechanism.  The  barrier  to  go  from  1  to  S  in  Fig.  18(a)  is  0.362 
eV,  but  the  barrier  from  S  to  1'  is  too  high  (2.557  eV). 

There  is  a  mechanism  for  transporting  a  particle  firom  the  side 
of  a  string  to  its  end,  and  it  turns  out  to  be  rather  surprising.  If  we 
push  the  atom  1  (Fig.  18(a))  towards  the  string,  1  will  replace  the 
atom  2  and  2  will  climb  on  the  string  and  form  a  trimer.  A  side  view 
of  this  trimer  is  shown  in  Fig.  18(b).  Once  placed  by  this  exchange 
mechanism  on  top  of  the  dimer  string,  atom  2  travels  rapidly  along 
the  string  and  descends  from  it  at  the  end,  as  shown  in  Fig.  19(a).  If 
the  string  has  an  F  site  vacancy  along  the  way,  the  atom  2  will  repair 
it  as  shown  in  Fig.  19(b).  Thus  the  exchange  mechanism  answers 
two  of  the  questions  posed  at  the  beginning  of  this  section:  it 
transports  the  atoms  from  the  side  of  the  string  to  its  end,  if  the  string 


The  exchange  process  takes  place  with  a  barrier  of  1.07  eV,  if  the 
atom  1  is  located  as  shown  in  Fig,  18(a),  or  with  a  barrier  of  1.36  eV  if 
atom  1  starts  from  the  site  S.  Both  these  barriers  are  rather  high  and 
the  exchange  is  not  very  efficient.  Fortunately,  the  exchange  rate  is 
increased  if  a  second  adatom  N  approaches  the  adatom  1  located  at 
the  side  of  the  string  (Fig.  20(a)).  They  form  an  1  -  N  dimer  whose 
configuration  is  shown  in  Fig.  20(b).  If  the  atom  N  in  Fig.  20(b)  is 
pushed  towards  the  string  the  atom  1  will  engage  in  the  exchange 
process;  the  outcome  is  shown  in  Fig.  20(c).  The  presence  of  N 
lowers  the  barrier  for  the  exchange  process  firom  1,07  eV  to  0.51  eV. 

The  exchange  process  is  not  the  only  option  available  to  the  1  -  N 
dimer  shown  in  Fig.  20(b).  It  can  also  move  away  firom  the  string  by 
breaking  the  weak  4-1  bond,  to  reach  the  configuration  shown  in  Fig, 
20(d).  The  barrier  for  this  move  is  0.34  eV.  The  1  -  N  dimer  will  be  a 
nucleation  center  for  the  growth  of  a  new  dimer  row  adjacent  to  the 
first.  Strings  consisting  of  double  and  triple  dimer  rows  have  been 
observed  experimentally.  Their  relative  frequency  depends  on  the 
growth  conditions. 

These  findings  provide  the  following  description  of  the  most 
likely  pathway  for  string  growth  and  error  repair.  Single  atoms 
moving  along  the  dimer  rows  of  the  substrate  have  no  trouble 
reaching  the  side  of  an  existing  string.  They  will  stick  to  it  and  are 
unlikely  to  jump  away  or  to  move  along  the  side.  There  is  a 
reasonable  probability  that  through  the  exchange  process  the  atom 
incorporates  itself  in  the  string,  and  lifts  one  of  the  string’s  atoms  on 
top  of  it.  Once  on  top  the  atom  moves  rapidly  along  the  string;  if  it 
encotmters  a  vacancy  it  falls  in  it  and  repairs  it;  if  not,  it  reaches  the 
end  of  the  string,  steps  down  from  it  and  attach  itself  to  the  string 
increasing  its  length.  Another  likely  scenario  is  that  the  atom  stuck 
t  the  side  of  a  string  might  wait  there  until  another  atom 
approaches  it.  The  presence  of  the  second  atom  aids  the  exchange 
process;  however  they  may  also  decide  to  move  jointly  away  from  the 
string  and  form  an  adjacent  dimer  string. 


Fig.  20.  (a)  An  adatom  at  the  side  of  a  dimer  row  before  engaging  in 
the  exchange  process,  (b)  The  position  N  of  a  second  atom,  which 
helps  1  engage  in  the  exchange  process,  (c)  The  atom  1  starting  in 
the  position  shown  in  (b)  undergoes  an  exchange  process,  and  put  the 
atom  4  on  top  of  the  dimer  row.  (d)  The  two  atoms  can  move  to  create 
a  nucleation  center  for  the  growth  of  a  second  dimer  row. 


6.  ADATOM  MOTION  NEAR  STEPS  AND  SINGLE  TO  DOUBLE 
STEP  TRANSITIONS 


Before  discussing  kinetics  on  stepped  Si(100)-2X1  terraces  we 
define  some  nomenclature  which  is  best  understood  by  looking  at  Fig. 
21.  This  shows  a  stepped  surface  with  one-atom-kigh  steps.  Such  a 
surface  has  two  kinds  of  terraces:  a  TA  terrace,  on  which  the  dimer 
rows  are  perpendicular  to  the  adjacent  steps,  and  a  TB  terrace  on 
which  the  rows  are  parallel.  The  TA  and  TB  terraces  alternate.  A 
step  having  a  TA  terrace  above  it  is  called  an  SA  step;  one  having  a 

TB  terrace  above  it  is  an  SB  step.^ 


Fig.  21.  A  vicinal  surface  with  one-atom-high  steps.  The  down  stairs 
direction  is  from  the  left  to  the  right.  The  arrows  show  the  fast 
pathways  of  adatom  migration  on  the  terraces. 


On  surfaces  having  two  atom  high  steps,  the  dimer  rows  on  all 
terraces  are  in  the  same  direction.  The  terraces  with  the  dimer  rows 
perpendicular  to  the  steps  (called  TB  terraces)  are  more  stable  than 
the  ones  having  dimer  rows  parallel  to  the  steps  (TA  terraces).  The 
steps  separating  the  TB  terraces  are  called  DB  steps.^ 

An  interesting  experiment  was  performed  by  Hoeven  et  al.,^ 
who  deposited  half  a  monolayer  of  silicon  on  a  Si(100)-2X1  surface 
having  S  steps  and  examined  the  result  by  STM.  They  found  that 
after  the  deposition  the  surface  had  DB  steps,  with  the  dimer  rows 
pexpendicular  to  the  steps  (  TB  terraces).  The  width  of  these  terraces 
is  twice  the  width  of  the  single  atom  high  terraces  of  the  initial 
surface. 

The  ability  of  the  adsorbed  atoms  to  create  such  a  highly 
organized  structure  by  undergoing  uncorrelated,  independent 
random  jumps  is  rather  striking.  In  this  section  we  look  for  the 
kinetic  pathway  by  which  this  transformation  takes  place.  As  a 
byproduct  we  also  find  a  possible  explanation  for  the  reason  why  the 
SB  steps  are  much  rougher  than  the  SA  ones.^^ 

The  overall  idea  is  rather  simple.  The  experiment  of  Hoeven  et 
al.  ^  starts  by  the  deposition  of  half  a  monolayer  of  Si  atoms  on  a 
surface  on  which  TA  and  TB  terraces  alternate  (Fig.  22(a)).  To  build 
the  double  layer  structure  shown  in  Fig.  22(b)  the  deposited  atoms 
must  migrate  in  a  peculiar  way:  the  ones  that  have  landed  on  a  TA 
terrace  must  stay  there;  the  ones  that  have  landed  on  a  TB  terrace 
must  move  onto  the  neighboring  TA  terraces  by  either  climbing  up  an 
SA  step  or  dovm  an  SB  step  (Fig.  22(c)).  When  half  a  monolayer  has 
been  deposited  the  TA  terraces  will  be  completely  covered.  The  length 
of  the  double  terraces  formed  in  this  way  is  tvdce  the  length  of  the 
initial  single  terraces. 

This  simple  model  makes  predictions  that  can  be  tested 
experimentally.  For  example,  if  a  quarter  of  a  monolayer  is  deposited 
the  surface  will  consist  of  single  atom  high  terraces;  the  length  of  the 
TB  terraces  will  be  roughly  equal  to  3L/2,  while  that  of  the  TA 
terraces  is  L/2.  Here  L  is  the  length  of  the  starting  terraces. 

We  proceed  now  to  find  if  the  mechanism  proposed  above  is 


substantiated  by  a  more  detailed  study.  To  do  this  we  determine  the 
binding  sites  at  the  step  and  at  the  terrace  edge,  compute  the  energy 
barriers  for  the  site  to  site  jumps  and  use  them  to  estimate  the 
jumping  rate  constants.  Then,  we  determine  how  the  deposited 
atoms  migrate  to  form  double  height  terraces  by  assuming  that  they 
execute  most  often  the  jumps  having  the  highest  rates. 


Fig.  22.  A  schematic  representation  of  the  kinetic  mechanism  for  the 
transformation  of  single-atom-high  steps  into  two-atom-high  steps. 


We  examine  now  the  steps  invoked  by  the  above  mechanism  for 
the  double  atom  step  formation.  First  we  look  at  the  motion  across 
the  TA  terrace.  According  to  the  above  mechanism  this  must  be  slow 
so  that  the  atoms  that  land  on  TA  will  have  a  high  probability  to  be 
stuck  there.  This  is  exactly  what  we  find.  We  have  already  shown  in 
the  previous  sections  that  an  adatom  will  jump  mainly  along  a  dimer 
row,  as  indicated  by  the  solid  arrows  in  Fig.  21.  Since  an  atom 
moving  along  the  dimer  row  is  likely  to  occupy  the  site  B,  E,  or  C,  we 
have  placed  an  atom  at  each  of  these  sites  respectively  and  moved  it 
towards  the  step.  The  fastest  pathway  we  foxmd  for  moving  across 
the  TA  terrace  encounters  a  barrier  of  1.25  eV.^®  This  is  higher  than 
other  relevant  barriers  and  most  of  the  atoms  landing  on  the  TA 
terrace  will  stay  on  it  and  form  islands,  just  as  assumed  in  the 
mechanism  proposed  above. 

We  must  now  see  if  we  can  find  a  rapid  pathway  for  the 
migration  of  an  atom  from  the  TB  terrace  across  the  bordering  steps 
and  onto  the  neighboring  TA  terraces.  We  discuss  first  the  motion 
across  the  SB  step  (Fig.  21). 

Since  on  a  TB  terrace  the  dimer  rows  are  perpendicular  to  the 
steps,  the  adsorbed  atom  has  no  difficulty  migrating  along  the  dimer 
row  and  reaching  the  SB  step.  The  only  possible  difficulty  may  be  the 
last  jump,  from  the  site  to  the  site  Bj^  on  the  dimer  nearest  to  the 
step  (  see  Fig.  21),  We  find  that  the  largest  barrier  encountered  by  an 
atom  moving  from  to  Bj^  is  0.44  eV.  This  is  substantially  larger 
than  the  0.29  eV  barrier  encoimtered  by  an  atom  moving  from  say 
to  B2  but  it  will  not  provide  'a  substantial  kinetic  handicap.  The 
terrace  is  very  narrow  and  an  atom  that  did  not  succeed  going  from 
Cj  to  Bj^  has  many  chances  to  tiy  ag£  .n. 

An  atom  that  has  reached  the  site  can  descend  across  the  SB 
step  onto  the  TA  terrace  and  reach  the  site  Ej  (  Fig.  21);  to  do  this  it 
must  overcome  two  barriers,  one  of  0.40  eV  and  the  other  of  0.84  eV. 
Once  at  Ej^,  it  can  no  longer  leave  the  step  edge:  the  barriers  for 
moving  back  across  the  SB  step  on  the  TB  terrace  or  away  from  E^  on 
the  TA  terrace,  are  both  of  about  2  eV.  This  result  is  consistent  with 
the  experimental  observation  of  Mo  and  Legally, 5  who  found  that  the 
SB  step  traps  the  adsorbed  particles  very  efficiently. 


On  the  other  hand,  the  atom  bound  at  can  easily  move  along 
tlie  SB  step.  The  largest  barrier  encountered  along  the  way  is  0.37 
eV.  This  is  expected  since  the  pathway  along  the  step  is  on  top  of  a 
dimer  row,  where  the  mobility  is  rather  high.  Thus,  an  atom 
descending  from  a  TB  terrace  across  the  SB  step  to  the  site  will 
move  along  the  step  until  it  encounters  another  adatom  on  the  same 
string  and  form  a  dimer.  This  dimer  will  line  up  with  the  dimer  row 
of  the  TB  terrace,  leading  thus  to  the  SB  step  roughening,  and  the  TB 
terrace  growth. 

The  descent  of  an  isolated  atom  from  a  B^  site  to  E]^  is  rather 
slow.  The  rate  of  this  transition  is  however  greatly  increased  by  the 
presence  of  another  atom  on  the  TA  terrace  near  the  step.  One 
example  is  shown  in  Fig.  23(a)  where  an  atom  at  E^^  helps  the  atom  at 
a  neighboring  Bj  site  to  jump  down  across  the  step.  The  two  atoms 
move  to  form  a  misaligned  dimer  (Fig.  23(b))  which  later  takes  the 
correct  position  (Fig.  23(c))  by  rotating  by  90  degrees.  The  highest 
barrier  in  this  sequence  is  0.73  eV. 

Next  we  examine  how  an  atom  on  a  TB  terrace  approaches  the 
SA  step  bordering  and  climbs  up  across  it  to  go  on  the  adjacent  TA 
terrace.  Such  an  atom,  denoted  1  in  Fig.  24(a),  reaches  the  position  G 
without  difficulty,  by  moving  on  top  of  a  dimer  row  as  indicated  by  the 
arrow.  Once  there,  the  atom  is  prevented  by  a  barrier  of  over  2  eV 
from  moving  back  onto  the  TB  terrace.  The  barrier  for  moving  along 
the  SA  step  to  the  S  site  is  0.36  eV.  Once  there  the  atom  is  trapped. 
Before  hopping  to  S  the  atom  has  a  chance  to  participate  in  an 
exchange  process.  The  outcome  of  this  process  is  shown  in  Fig.  24(b): 
the  atom  1  replaces  the  atom  2,  and  the  latter  is  lifted  up  on  the  TA 
terrace.  This  exchange  process  must  overcome  an  energy  barrier  of 
1.07  eV.  The  rate  is  increased  greatly  if  another  adatom,  denoted  N 
in  Fig.  24(c),  is  located  nearby.  The  precursor  configuration  for  the 
exchange  is  shown  in  Fig.  24(d).  If  we  push  the  atom  N  towards  the 
step  the  atom  1  will  engage  in  an  exchange  process  whose  outcome  is 
shown  in  Fig.  24(e).  The  presence  of  N  lowers  the  barrier  for  the 
exchange  process  from  1.07  eV  to  0.51  eV. 

The  atom  placed  by  the  exchange  process  on  top  of  the  dimer  row 
(on  the  TA  terrace)  can  easily  move  in  a  direction  parallel  to  the  step. 


leaving  behind  the  atom  N  in  the  configuration  shown  in  Fig.  24(a). 


I 


Pig.  23.  The  manner  in  which  the  migration  of  the  atom  1  across  the 
SB  step  is  aided  by  the  presence  of  the  atom  2  near  the  step  (on  the 
lower  terrace). 

The  exchange  process  is  not  the  only  possible  move  taking  place 
from  the  configuration  shown  in  Fig.  24(d):  the  atom  N  can  break  the 
weak  1-N  bond  and  move  away  onto  the  TB  terrace.  The  barrier  for 
this  is  0.9  eV. 

Overall,  the  SA  step  is  less  effective  than  the  SB  step  in  trapping 


the  adatoms  reaching  it  from  the  lower  terrace.  This  is  consistent 
with  the  observation  of  Mo  and  Legally^®  that  more  islands  are 
formed  on  a  TA  terrace  near  the  SA  step  than  on  the  TB  terrace  near 
the  SA  step. 
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Fig.  24.  The  exchange  mechanism  by  which  an  atom  deposited  on 
the  TB  terrace  is  transferred  onto  the  SA  terrace  by  climbing  up  the 


SA  step. 

We  can  now  summarize  the  kinetic  factors  that  lead  to  the 
transformation  of  single-atom-high  steps  into  double-atom-high  ones 
after  the  deposition  of  half  a  monolayer  of  silicon  on  the  surface.  The 
first  is  the  extreme  anisotropy  of  adsorbate  motion  on  the  Si(100)-2X1 
surface.  An  adatom  located  on  top  of  a  dimer  row  moves  frequently 
along  the  row,  and  rarely  in  the  direction  perpendicular  to  it.  This 
affects  strongly  the  particle  transport  on  terraces.  An  atom  adsorbed 
on  a  TB  terrace,  on  which  the  dimer  rows  are  perpendicular  to  the 
step,  has  no  trouble  reaching  a  terrace  edge.  On  a  TA  terrace  the 
dimer  rows  run  along  the  terrace  and  the  transport  of  the  adsorbed 
particles  to  the  edge  or  the  step  is  slow.  This  effect  alone  will  make  it 
difficult  for  the  adsorbed  atoms  to  leave  the  TA  terraces. 

The  einisotropy  of  diffusion  is  only  half  the  story.  It  explains  the 
tendency  of  the  adsorbed  particles  that  have  landed  on  a  TA  terrace  to 
stay  on  the  same  terrace.  To  understand  the  formation  of  the  double 
layer  steps  we  must  also  find  why  and  how  the  particles  landing  on 
the  TB  terraces  leave  them  and  move  to  the  TA  ones.  Our 
calctdations  indicate  that  a  particle  can  leave  a  TB  terrace  by 
climbing  down  (across  the  SB  step)  onto  the  neighboring  TA  terrace; 
the  largest  barrier  along  the  way  is  0.84  eV.  Furthermore,  the  atoms 
adsorbed  on  the  TB  terraces  can  reach  the  SA  step  and  engage  in 
what  we  call  the  exchange  process;  as  a  result  an  adatom 
disappearing  from  the  TB  terrace  reappears  on  the  TA  one.  The 
barrier  for  the  exchange  process  is  1.07  eV  if  imassisted  by  a  second 
atom,  and  0.54  eV  if  a  second  atom  is  nearby. 


7.  KINK-KINK  INTERACTIONS  IN  STEP  ROUGHENING 

It  has  been  observed  experimentally  that  the  single-atom-high  B 
steps  (SB)  on  a  Si(100)-2X1  surface  has  kinks,  while  the  SA  steps  are 
practically  straight.^^  In  Sec.  6  we  have  proposed  that  the  step  is 
rough  because  to  a  good  approximation  the  step  is  formed  by  the 
independent  growth  of  the  dimer  rows  of  the  TB  terrace.  Here  we 
present  a  thermodynamic  study  of  roughening  based  on  a 
phenomenological  model  for  the  energy  of  the  kinks  and  the 


interaction  between  them.^®  Understanding  and  controlling  the 
shape  of  the  steps  is  especially  important  because  it  may  affect  the 
quality  of  the  tilted  superlattices  which  are  grown  on  stepped 
surfaces. 

The  work  that  is  most  significant  to  this  section  is  that  of 
Swartzentruber  et  who  studied  the  statistical  properties  of  the 

rough  SB  steps  in  the  STM  measurements.  They  noted  that  these  can 
be  reproduced  rather  well  by  a  model  that  assiimes  that  the  energy  to 
form  a  step  is  of  the  form  a  +  b  n,  where  n  is  the  number  of  atoms  in  a 
kink  (we  give  below  a  more  precise  definition)  and  a  is  a  corner 
formation  energy.  Here  we  propose  a  slight  extension  of  this  work: 
we  show  that  an  energy  expression  that  includes  a  kink-kink 
interaction  term  gives  a  slightly  better  fit  of  the  data.  Unfortunately 
there  is  insufficient  data  to  decide  between  these  two  models.  Future 
measurements  of  the  statistical  properties  of  the  kinks  as  a  function 
of  temperature  will  be  very  helpful  in  determining  the  magnitude  of 
the  kink-kink  interactions. 

The  properties  of  the  kinks  are  characterized  by  two  distribution 
functions:  the  probability  N(n)  of  finding  a  kink  of  height  n,  and  the 
probability  P(s)  that  the  distance  between  two  neighboring  kinks  is  s. 
The  variables  s  and  n  are  defined  in  the  insert  of  Fig.  25(a). 
Swartzentruber  et  have  determined  these  two  functions  fi:om 

the  STM  pictures  of  the  B  and  A  steps  on  a  Si(100)-(2X1)  sxirface.  They 
interpreted  the  results  by  using  an  independent  kink  (IK)  model 
which  gives 

N(n)  -  expC-E(n)/kT],  (7.1) 

P(s)  =  p(l-p)®^-^.  (7.2) 

Here  T  is  the  absolute  temperature,  p  *  0.44  is  the  probability  that  a 
kink  exists  at  a  site  along  the  step,  and  E(n)  =  e  n  +  C,  where  C  is  the 
energy  of  the  two  corners  in  a  kink.  Fitting  the  data  for  a  SB  step 
gives  e  =  0.028  ±  0.002  eV  and  C  =  0.08  ±  0.02  eV. 

The  IK  model  fits  the  data  fairly  well.  There  are  however  small 
discrepancies:  (a)  at  low  values  of  s  the  measured  P(s)  deviates 
considerably  from  Eq.  (2);  (b)  for  large  n  (  ^  16)  the  measured  values 
of  N(n)  are  higher  than  the  values  given  by  the  IK  model;  (c)  in  the 


strictly  independent  kink  model  if  one  walks  along  the  step  one  would 
have  equal  probability  to  wonder  upwards  or  downwards,  while 
experimentally  the  directions  of  the  roughening  at  two  neighboring 
sites  are  correlated;  (d)  the  amplitude  of  the  roughening  in  a  step 
generated  by  using  the  IK  model  with  the  parameters  given  above,  is 
smaller  than  that  of  the  experimental  pattern. 
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Fig.  26.  (a)  The  kink  height  distribution,  (b)  The  kink-kink  distance 
distribution.  The  data  from  Ref.  22  is  represented  by  stars,  the 
soliddots  and  the  triangles  give  the  distributions  calculated  with  and 
without  kink-kink  interactions,  respectively. 


Several  recent  articles  have  used  phenomenological 
Hamiltonians  to  study  the  equilibritam  structures  of  Si(lOO)  stepped 
surfaces.21.62-70  in  most  of  these  studies  the  interaction  between 
steps  has  been  taken  into  accoimt,  while  the  effects  of  the  kink-kink 
interaction  along  a  step  have  not  been  treated  explicitly.  Since  the 
kink-kink  distance  can  be  much  shorter  along  the  same  step  than 
step-step  separation,  the  interaction  between  the  kinks  may  have 
stronger  effects. 

In  what  follows  we  use  the  Hamiltonian 

H=II  [(Injl  e  +  C)5ij  -D(l-  5ij)ni  nj  /  (7.3) 

Hi  nj 

which  adds  to  the  IK  expression  a  kink-kink  interaction.  The  kink 
height  nj  and  the  kink-kink  distance  s^j  are  defined  in  Fig.  25(a).  The 
kink-kink  coupling  strength  D  is  positive.  The  sign  of  the  product 
UjUj  is  given  by  the  following  rule:  walk  along  the  step  starting  from 
one  end  and  mark  with  arrows  the  direction  in  which  you  walk  along 
the  kink;  each  kink  is  thus  assigned  an  arrow  and  if  the  arrows  for 
the  kinks  i  and  j  are  antiparallel  (see  Fig.  26(a))  njnj<0  ;  if  they  are 
parallel  (  see  Fig.  26(b))  n|nj>0.  The  kink-kink  interaction  gives 
preference  to  configurations  in  which  the  arrows  of  the  neighboring 
kinks  are  parallel;  thus  it  favors  deeper  kinks  and  smaller  kink-kink 
separations. 


I 
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Fig.  26. 


We  use  the  (sy)’^  dependence  because  it  is  the  one  that  gives  the 
best  agreement  with  the  experimental  data.  We  have  also  tried  the 
forms  n^ajS^j  and  s^j  and  found  them  less  successful. 

There  are  some  arguments  in  the  literature^  ^  that  favor  a  (sy)'^ 
dependence,  but  we  do  not  find  them  compelling  and  prefer  to  regard 
our  Hamiltonian  as  empirical.  A  convincing  derivation  of  the 
effective  Hamiltonian  would  obviously  be  valuable. 

For  njnj>0  the  Hamiltonian  Eq.  (7.3)  favors  kinks  with  infinite 
height.  In  reality,  the  maximum  kink  height  is  limited  by  the  terrace 
width  L,  so  this  quantity  should  in  principle  be  a  parameter  in  the 
theory.  While  adding  a  term  limiting  the  maximum  kink  height, 
such  as  the  one  expressing  step-step  interaction,  would  not  greatly 
complicate  the  simulations  below,  we  find  that  we  can  fit  well  the 
existing  data  with  the  Hamiltonian  (7.3).  This  happens  because,  at 
the  temperature  and  angle  of  cut  at  which  the  data  was  taken,  the 
first  term  in  Eq.  (7.3)  discriminates  effectively  against  the  creation  of 
kinks  with  high  values  of  n. 

The  effects  of  the  interaction  term  can  be  examined  by  first 
generating  kinks  on  a  step  by  Monte  Carlo  (MC)  simulations  and 
then  calculating  the  distribution  functions  N(n)  and  P(s).  We  start 
with  an  ideal  step  of  length  300  X  (2a)  where  2a  is  the  period  of 
translational  symmetry  in  the  direction  parallel  to  the  step,  and  a  is 
the  surface  lattice  constant  on  the  Si(lOO)  surface  before 
reconstruction.  Periodic  baundary  conditions  in  the  direction 
parallel  to  the  step  are  used  with  the  period  600  a.  The  MC  program 
chooses  randomly  a  site  along  the  step  and  attempts  to  remove  two 
adjacent  dimers  from  the  step  and  place  them  at  another  randomly 
chosen  site  along  the  step.  The  Hamiltonian  Eq.  (7.3)  is  used  to 
determine  whether  such  a  change  is  accepted.  It  takes  about  200,000 
moves  to  reach  the  equilibrium  kink  configuration.  To  test  whether  a 
high  n  cut  off  is  necessary,  to  prevent  the  formation  of  kinks  whose 
height  exceeds  the  width  of  the  terrace,  we  have  rejected  in  our 
simulations  attempts  to  create  kinks  with  n>25  X  (2a).  The  munber  of 
attempts  rejected  by  this  restriction  is  negligible,  which  indicates  that 
the  tendency  of  the  kink-kink  interaction  to  favor  high  n  values  is 


statistically  insignificant. 

We  move  two  dimers  each  time  as  the  removal  of  two  dimers  is 
energetically  favored:  it  allows  the  dimerization  of  the  atoms  left 

behind.  34 

In  our  simulations  we  have  systematically  examined  the  effects 
of  the  interaction  term  by  var3dng  the  magnitude  of  D.  Here  we  only 
compare  the  best  fits  of  the  distributions  N(n)  and  P(s)  given  by  the 
present  model  and  the  IK  model,  at  T=600  ®C.  These  are  shown  in 
Fig.  25.  The  best  fitting  parameters  are  e  =  0.029  eV,  C  =  0.050  eV, 
and  D=0.0052  eV,  for  the  present  model,  and  e  =  0.028  eV,  C  =  0.063 
eV,  and  D=0.0,  for  the  IK  model.  The  corner  energy  C  =  0.063  eV 
obtained  by  us  is  slightly  lower  than  the  value  recommended  by 
Swartzentruber  et  The  present  model  agrees  better  with  the 

measurements,  especially  for  large  values  of  n  and  small  values  of  s. 
We  have  also  calculated  these  distributions  for  T=500  °C,  by  assuming 
that  the  constants  in  the  Hamiltonian  are  temperature  independent, 
and  found  that  they  do  not  change  significantly. 

The  kink-kink  interaction  has  a  greater  influence  on  the  kink- 
kink  distance  distribution  P(s)  for  smaller  s  (Fig.  25(b)).  The 
parameter  set  {e,  C,  D)  that  improves  the  fitting  of  N(n)  also  gives  a 
better  fitting  of  P(s).  In  particular,  the  inclusion  of  the  interaction 
term  increases  P(2)  as  compared  to  the  value  given  by  the  IK  model 
and  diminishes  the  magnitude  of  P(s)  for  s  ^  4.  These  changes 
improve  the  agreement  with  the  experiment. 

In  Fig.  27,  we  also  compare  the  step  patterns  obtained  for 
different  values  of  D,  for  a  terrace  width  L=34  X  (2a),  corresuonding  to 
an  angle  of  cut  of  0.3°  towards  [100].  As  expected,  the  inclusion  of  the 
interaction  increases  thermal  roughening,  and  also  influences  the 
global  pattern  of  the  step  to  favor  n^nj  >  0.  The  SB  step  structure 
shown  in  Fig.  27(a)  (for  D  =  0.0052  eV)  resembles  more  closely  Fig.  1  of 
Ref.  22,  than  the  pattern  for  0=0  (Fig.  27(b)),  especially  when  the 
magnitude  of  the  thermal  roughening  is  considered. 
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